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Abstract 

We present an analytic five-loop calculation for the additive renormaliza- 
tion constant A{u, e) and the associated renormalization-group function B(u) 
of the specific heat of the 0(n) symmetric (/>^ theory within the minimal 
subtraction scheme. We show that this calculation does not require new five- 
loop integrations but can be performed on the basis of the previous five-loop 
calculation of the four-point vertex function combined with an appropriate 
identification of symmetry factors of vacuum diagrams. We also determine 
the amplitude function of the specific heat in three dimensions for 

n = 1,2,3 above Tc and F^(u) for n = 1 below Tc up to five-loop order, 
without using the e = 4 — d expansion. Accurate results are obtained from 
Borel resummations of B{u) for n = 1,2,3 and of the amplitude functions for 
n = 1. Previous conjectures regarding the smallness of the resummed higher- 
order contributions are confirmed. Combining our results for B{u) and F^(u) 
for n = 1, 2, 3 with those of a recent three- loop calculation of F^{u) for gen- 
eral n in d = 3 dimensions we calculate Borel resummed universal amplitude 
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ratios A'^jA' for n = 1, 2, 3. Our result for A^^jA' = 1.056 ±0.004 forn = 2 
is significantly more accurate than the previous result obtained from the e 
expansion up to O(e^) and agrees well with the high-precision experimental 
result A^/A~ = 1.054 it 0.001 for ^He near the superfluid transition obtained 
from a recent experiment in space. 
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I. INTRODUCTION 



One of the fundamental achievements of the renormahzation-group (RG) theory of crit- 
ical phenomena is the identification of universality classes in terms of the dimensionality d 
of the system and the number n of components of the order parameter . Specifically, RG 
theory predicts that the critical exponents, certain amplitude ratios and scaling functions 
are universal quantities that do not depend, e.g., on the strength of the interaction or on 
thermodynamic variables (such as the pressure). The superfluid transition of ^He belongs 
to the d = ?>, n = 2 universality class and provides a unique opportunity for an experimen- 
tal test of the universality prediction by means of measurements of the critical behavior at 
various pressures P along the A-line T\{P). Early tests have been performed by Ahlers and 
collaborators and consistency with the universality prediction was found within the experi- 
mental resolution 0. At a significantly higher level of accuracy, the superfluid density and 
the specific heat (or, equivalently, thermal expansion coefficient) above and below Tx{P) are 
planned to be measured in the Superfluid Universality Experiment (SUE) P] under micro- 
gravity conditions or at reduced gravity in the low-gravity simulator ||^. As demonstrated 
recently 0, this would allow to perform measurements up to |t| ~ 10^^ in the reduced 
temperature t = {T - Tx{P)) /Tx{P). 

On the theoretical side, the corresponding challenge is to calculate as accurately as 
possible the properties of the 0{n) symmetric 0^ model in three dimensions. To extract the 
leading critical exponents from the experimental data and to demonstrate their universality 
at a highly quantitative level requires detailed knowledge on the ingredients of a nonlinear 
RG analysis 0. They include not only the well-known RG exponent functions of the 0^ 
model whose fixed point values determine the critical exponents but also the less well-known 
amplitude functions [[^-[T2[| which contain the information about universal ratios of leading 
and subleading amplitudes 

The existing theoretical predictions on the critical exponents within the minimal 



subtraction scheme p^,p!5| are based on field-theoretic calculations to five-loop order [|T^-[T9 
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and Borel resummation. By contrast, the present theoretical knowledge of the amplitude 
ratios for n > 1 below is based only on low-order (mainly 1- and 2-loop) calculations 
which imply an uncertainty at the level of at least 10-30% [|1|. It has therefore been proposed 
pO| to significantly reduce this uncertainty by performing new higher-order field-theoretic 
calculations and Borel resummations of various amplitude functions in three dimensions. 

Both conceptual and computational steps towards this goal have already been performed. 
The conceptual progress includes the demonstration that the d = 3 field theory suggested 
by Parisi |2ll] can well be realized within the minimal subtraction scheme at d = 3 



by incorporating Symanzik's non-vanishing mass shift p2[ and that spurious Goldstone 
singularities for n > 1 below can well be treated within this approach [0 by using an 
appropriately defined pseudo-correlation length 0. The computational steps include the 
determination of the amplitude functions F^{u) and F^{u) of the specific heat in three 
dimensions above Tc for n = 1,2,3 [|T0| and below Tc for n = 1 [|ll|, respectively, up to 



five-loop order, and their Borel resummation. These calculations, however, were not yet 
complete because of an approximation regarding the additive renormalization A{u, e) of the 
specific heat and the associated RG function B{u). Due to the lack of knowledge in the 
literature about higher-order terms, A{u,e) and B{u) were approximated by their two-loop 
expressions. Although the good agreement between low-order d = 3 perturbation results 
||7| j23|j2^ and accurate experiments |0,|25|j2^ provided some indication for the smallness of 



the effect of the higher-order terms of B{u), no reliable estimate could be given for the 
remaining uncertainty of -F±(m) which could well be of relevance at the level of accuracy 
anticipated in future experiments 0. Furthermore we recall that any inaccuracy of B{u) 
enters not only the formulas ||^ for several universal amplitude ratios but also the formulas 
needed to determine the effective coupling u{l) from the specific heat |7| J23|J2^ . 



It is the purpose of the present paper to provide the missing information on the higher- 
order terms of A{u,6) and B{u) by means of a new five-loop calculation. We shall show 
that the analytic calculation of A{u,e) and B{u) can be directly related to the previous 
calculations [[TBHTSl] of the four-point vertex function. This provides the crucial simplification 



that no new evaluations of three-, four- and five-loop integrals are necessary but that only 
a new determination of symmetry and 0{n) group factors of vacuum diagrams is sufficient. 

Using the five- loop expression of B{u) we are in the position to determine the correct 
higher-order terms of the minimally renormalized amplitude functions for n = 1, 2, 3 

and F_{u) for n = 1 in three dimensions on the basis of previous work where a 

different renormalization scheme was used. The new coefficients of the higher-order terms of 



-F+(m) turn out to differ considerably from the previous approximate coefficients [|TOl whereas 
the coefficients of F^{u) are only weakly affected by the new higher-order terms of B{u). 

We also perform new Borel resummations of B{u) for = 1,2,3 as well as of F-{u) 
and of F_{u) — F^{u) for n = 1. It turns out that the result of the Borel resummation for 
B{u) including the new terms up to five- loop order differs from the two- loop result B{u) = 
n/2 + 0{u^) by only less than 1% at the fixed point. For the amplitude functions, our new 
Borel resummation results differ from the previous ones [p!0| , p!I| by about 1% for F_ and by 
less than 0.1% for F_ — F^ at the fixed point. This is a nontrivial and important confirmation 
of the previous conjectures about the smallness of resummed higher-order contributions 

Ml- 

As a first application, we calculate the universal ratios A'^/A" and /a~ of the leading 
and subleading amplitudes of the specific heat above and below for n = 1 and d = ?>. In 
addition we calculate Borel resummation values of A'^ /A^ for n = 2 and 3 by combining 
our present results for B{u) and F^{u) for n = 2 and 3 with those of a recent three- loop 



calculation of F_{u) for general n ^T]]. All of our calculations are performed at = 3 
dimensions, without using the e expansion. Our result A'^ /A'^ = 1.056 ± 0.004 for n = 2 
is more accurate than the previous result 1.0294 ± 0.0134 obtained from the e = 4 — d 
expansion up to O(e^) and agrees well with the high-precision experimental result ^ 1.054± 
0.001 for ^He near the superfiuid transition obtained from a recent experiment in space. 
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II. ADDITIVE RENORMALIZATION OF THE SPECIFIC HEAT 



The 0(n) symmetric 0'* model is defined by the usual Landau-Ginzburg- Wilson func- 
tional 

7i{0o(x)} = d'x Qro0^ + ^ E(Wo.)' + uoi^y - h ■ 0o) (2.1) 

for the n-component field 0o(x) = (0oi(x), . . . , 0on(x)) where 

ro = roc + aot, t = {T-T,)/T,, (2.2) 

and hf) = {ho, 0, . . . , 0). The Gibbs free energy per unit volume (divided by ksT) is 

Fo(ro, uo, ho) = In J P0o expi-H) . (2.3) 

We shall consider the bulk limit V ^ oo. We are interested in the specific heat per unit 
volume at vanishing external field = (divided by Boltzmann's constant ks) where ± 
refers to T > Tc and T <Tc, respectively. Near Tc, is determined by 

C^ = Cb- T^^Foiro, Uo, 0) = Cb - ao^i^o(ro, Uo, 0) (2.4) 

where is an analytic background term. Alternatively the Helmholtz free energy per unit 
volume ro(ro, uq, Mq) = F{ro, uq, ho) + h^Mo with Mq = {<t>m) determines in the /iq — > 
limit according to 
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d 

= Cb~ ao^r^ro (ro, mq, Mo(ro, uo)) • (2.5) 

The perturbative expression for ro(ro, Uq, Mq) is obtained from the negative sum of all one- 
particle irreducible (1 PI) vacuum diagrams |]TB[. The perturbative expression for is 
then determined by the vertex functions f|^'°'' = d^Fo/drQ which we consider as functions 
of appropriately defined correlation lengths ^+ and above and below Tc @JT3, 



C^ = CB-altt^'\i±,uo,d). (2.6) 
A description of the critical behavior requires to turn to the renormalized vertex functions 



72f(2,0) 



(2.7) 



We work at infinite cutoff using the prescriptions of dimensional regularization and minimal 
subtraction at fixed dimension 2 < d < 4 without employing the e = 4 — d expansion 
The Z-factors are introduced as 



r = ^(ro - roc), u = fi ^A^Z^ ^Z^uq, 



(2.8) 



where the geometric factor 



Ad = r(3 - d/2) 2^-'^7T-'^/\d - 2)-i 



(2.9) 



becomes ^3 = [An) ^ for d = 3 and A4 = (Stt^) ^ for d = 4. These Z-factors Zi{u,e) and 
the associated field-theoretic functions M 



Criu) = n (9^1nZ^(M,e)"^ ^ , 
C4>{u) = nd^\nZ^{u,e)-^\^ , 
(3u{u, e) = -eu + /3{u) =u -e + /x df,{Z~'^Zl] 



(2.10) 
(2.11) 
(2.12) 



are known up to five- loop order [17-19 . 



The main quantity of interest in the present paper is the renormalization constant A{u, e) 
in Eq. ( p.TD which absorbs the additive poles of both f and f Previously |10|JTT[] 
A{u, e) was employed only in its two-loop form ^ 

1 u 

A(u, e) = -2n- - 8n(n + 2)- + 0(u^) . (2.13) 
e e"' 

Here we report on a calculation of A{u, e) up to five-loop order. We would like to stress that 
this calculation does not require new five-loop integrations but can be performed on the 
basis of the previous five-loop calculation []T7|-|T^ of the four-point vertex function combined 
with an appropriate identification of symmetry and 0{n) group factors of vacuum diagrams 
which are shown in Fig. |I]. Their negative sum determines the Helmholtz free energy Fq 
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up to five-loop order. For the present purpose of determining the pole terms at d = 4 it 
suffices to consider only the case tq > and Mq = where only four-point vertices exist. 
The diagrams are labeled (1) in one-loop order, (2) in two-loop order, (3), (4) in three- loop 
order, (5)-(8) in four-loop order and (9)-(18) in five-loop order. The analytic expression 
of an m-loop diagram (i) is given by the product of the coupling (— mq)™"^, the symmetry 
factor S^^\ the 0{n) group factor G^^\n) and the momentum integral expression I^^\rQ,e). 
Thus the structure of the diagrammatic expression of a typical diagram, e.g., (16) is 



5(i6)(-Mo)'G(^'^)j(i^)(ro,e) (2.14) 
2592 {-uoY 




4 n"^ + 8n^ + 32^2 + 40n 



81 

X / / / / / Gi ■ G2 ■ Gi+2-3 ■ G3 ■ Gi+2-4 ■ Gi ■ 6*1+2-5 ■ (2.15) 

with /p = {2tt)~'^ J d'^p and the propagators Gi±j = {tq + |pj ± Pjp)~^. The symmetry and 
group factors are listed in Table |. 

To calculate the additive renormalization constant A{u, e) one needs to calculate those 
ultraviolet d = 4 pole terms of the diagrams contributing to T^q^^ that are left after sub- 
traction of subdivergences. One obtains Fq^^"* by taking two derivatives of Fq with respect 
to tq. The analytic calculation of the poles of the diagrams for Fq^^'' is identical to that car- 
ried out previously [|T^-[T^ for the four-point vertex function tI^'^^K To see this, one should 



take into account that in the minimal subtraction scheme |T^ the ultraviolet pole terms 



specified above do not depend on tq. Then by using the method of infrared rearrangement 
||33| one can nullify ro and introduce for each diagram a new fictitious external momentum 
to regularize infrared divergences. Then one can see that only a particular subset of those 
diagrams of Fq"'^'' are relevant in the present context, namely those where the four external 
legs are connected to each diagram through only two four-point vertices (rather than three 
four-point vertices or four four-point vertices). The details of the calculation are presented 
in Appendix A. 
The result reads 
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(2.16) 



m=l 



where A^^^ denotes the contribution of m-loop order, 



A^'\u,e) 



3 _ 40 24(n + 4) 



e 



--n{n + 2) 
— n{n + 2) 



248n + 1024 A8{n + A){n + 5) 



u 



(2.17) 



(n + 8)(12C(3) -25) 96n + 696 



u 



(2.18) 
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?2(?2 + 2) 



768(n + 4)(n + 5)(5n + 28) 128(293^2 + 2624n + 5840) 



+ 



9216C(3)(5n + 22) + 32(519^^ + 8462n + 25048) 

192C(3)(7n2 - 28n + 48) + 4608C(4)(5n + 22) + 64(31^2 + 2354n + 9306) 



+ (48C(3)(3n2 - 382n - 1700) + 288C(4)(4n2 + 39n + 146) 
-3072C(5)(5n + 22) - 3(319^^ - 13968n - 64864)' ^ 



u 



(2.19) 



where C(s) = Ej°^iJ~' is the Riemann zeta function with C(3) = 1.20205690, C(4) = 7rV90 
and C(5) = 1.03692776. Most important is the (i-independent RG function B{u) which is 
determined by 

dA{u, e) 



AB{u) = [2Cr{u) - e] A{u, e) + (3u{u, e)- 



du 



(2.20) 



Using A{u, e) of Eqs. (p.l6|) -( p.l9|) and the perturbative expressions for (r and of Refs. 
l3-|19l we find 



n 



B{u) = - + 3n(n + 2)u^ - -n{n + 2)(n + 8)(25 - 12C(3))m'^ 
2 3 

+ ]^n{n + 2) [l6C(3)(3n2 - 382r2 - 1700) - 1024C(5)(5n + 22) 



+96C(4)(4n2 + 39n + 146) - 319n2 + 13968n + 64864 



u 



(2.21) 



The terms of 0{v?) and 0{v?') agree with those of Ref. [Q. In Table |I| the coefficients 
of the power series 
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oo 



B{u) = E c,„u" (2.22) 

m=0 

are given for n = 0,1,2,3 up to m = 4 corresponding to five- loop order. Table y also 
contains the coefficients /f™''' of the power series of the functions 

oo 

Mu) = E (2-23) 

m=l 

where fi{u) denotes the functions P{u), Cr{u) and ((f,{u) for i = 1,2,3, respectively. These 
coefficients are taken from Refs. pT|-[TP[|. Up to four-loop order they agree with those in 
Table 1 of Ref. §. (Note that ^ in the table caption of Ref. § should read /f ^ x 10"^) 
The five-loop coefficients fl^\ f^P and f^^ differ from those in Table 1 of Ref. according 
to the corrections in five-loop order in Ref. [|l^. 

In Fig. H the partial sums of B{u) from two- to five- loop order are shown for the example 
n = 2. As expected, the contributions for m > 2 have alternating signs and increase 
considerably in magnitude. Clearly a resummation of B{u) is necessary similar to that for 
(r{u), ((i,{u) and /3(m) performed previously [§. 

First we reexamine the fixed point values u*, (3u{u* , 1) = 0, for n = 1, 2, 3 obtained in 
Refs. P,|2^ by means of Borel resummation on the basis of previous five-loop results p7|Jl8[] 
and in Ref. [0 on the basis of four-loop results. Here we employ the corrected five-loop 
coefficients for the e expansion of the fixed point value which we have derived from Eq. (8) 
of Ref. [0. Employing the standard Borel resummation method [^p5| we have obtained 



the fixed point values in three dimensions 



= 0.0404 ± 0.0003 for n = 1 , (2.24) 
u'' = 0.0362 ±0.0002 for n = 2 , (2.25) 
= 0.0327 ± 0.0001 for n = 3 . (2.26) 



The corresponding resummation parameters a and h = 5.5 + n/2 |^ are 

2.22 <a< 3.41 , 6 = 6.0 for n = 1 , (2.27) 
2.45<a<3.43, 6 = 6.5 for n = 2 , (2.28) 
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2.71 < a< 3.43 



7.0 



for n = 3 . 



(2.29) 



The previous fixed point values [^,0,^] are consistent with Eqs. ( |2.24| )-( |2.26| ) within the 
previous error bars. The present error bars are smaller than the previous ones P, p!6| , pl . 
(The range of a determines our error bars, as described further below.) 

We have performed Borel resummations of B{u) at the fixed point u* for the cases n = 1, 
2, 3. In addition, for the important case n = 2 (superfiuid ^He), we have determined the 
Borel resummed function B{u) at various values of u. The results are given in Eq. (p.34|) - 
( p.36| ) and in Figs. || and |^. 



A description of the Borel resummation method for the present purpose has been 
given in Section 5 of Ref. 0. (In Eq. (5.10) of Ref. ajk should read aj^mM-) In the 
present work, however, we use a different way of determining the parameters a and b of the 
summations. This implies a different determination of the error bars. 

For B{u) the value of the parameter b is not known from an analysis of the large-order 
behavior (see Eq. (5.6) of Ref. [§] and references therein). Here we fix both b and a by 
requiring fastest convergence of the series of the partial Borel sums S"*-^^ = S^^\u,a,b) for 
B{u) defined in Eq. (5.12) of Ref. |^ (here L corresponds to {L + l)-loop order). To do so 
we look for the minima of A*^^^ and A^^-* with regard to variations of both b and a where 



A^^\u,a,b) 



(2.30) 



This yields five-loop values of the parameters a, b for each u. In order to define an error bar, 
we apply the same method to the four- loop result of B{u), i.e. to A*^^) and A*^^\ The four- 
loop values of a, 6 together with the five-loop values provide the ranges of the best values of 
a and 6 as a result of the combined four- and five-loop analysis. Then we define the error 
bar of the five- loop result for B{u) by the maximum and minimum of the resummed four- 
and five-loop values for B{u) over the ranges of the best values of a, b. At the fixed point 
u*, Eqs. ( |2.24D -( P?^ , we find the ranges 



0.95<a<1.08, 5.7<6<7.75 forn = 1 



(2.31) 
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0.94<a<1.04, 7.0<6<8.59 forn = 2 , (2.32) 
0.94<a<1.02, 8.18 <6< 9.76 for n = 3 . (2.33) 

The corresponding Borel resummed results for B{u*) are 

B{u*) = 0.5024 ± 0.0011 for n = l, (2.34) 
B{u*) = 1.0053 ± 0.0022 for n = 2 , (2.35) 
B{u*) = 1.5080 ± 0.0034 for n = 3 . (2.36) 

We have also determined the function B{u) for n = 2 (superfluid *^He) in the range < m < 
0.04 as shown in Fig. ^. 

Most remarkable is the smallness of the deviation of the resummed function B{u) for 
n = 1,2,3 from its two-loop approximation n/2. This confirms previous conjectures [P-[TT[| 
and justifies earlier analyses 0,^,^ . 



III. AMPLITUDE FUNCTIONS F± IN THREE DIMENSIONS 

A. Definition of F± 

By means of the renormalized vertex functions in Eq. ( ^.Tj ) we define the dimensionless 
amplitude functions 

F±{fi^±,u,d) = -4/iM^^ri''°)(e±,M,/i,ci). (3.1) 

They enter the critical behavior of the specific heat in three dimensions in the form of the 
functions 

F±(1,m,3) = F±(m) (3.2) 

according to 

C^ = Cs + ]a'f,-'AsK^{u{l^)) exp f^'"^ ^TT^cln' (3.3) 

4 Ju Pu[U , L) 
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where 

K^{u) = F^{u)-Aiu,l) (3.4) 

and a = Zr(u, l)~^ao. In Eq. (|0|), u{l) is the effective couphng satisfying 

l^ = Pu{uil),l) (3.5) 

with u{l) = u. The flow parameters /+ and /_ are chosen as /+ = and /_ = (/i^_)~^ 

above and below Tc. 

B. Power series of F± 

The amplitude functions are expandable in integer powers of m 0] and have the power 
series [fOU 



oo 



F+iu) = Y: otrX" (3.6) 



m=0 

and ITT 



-| oo 
" m=0 

We have determined cj„ up to five-loop order (i.e., up to m = 4) for n = 1, 2, 3 and 
up to five- loop order (i.e., up to m = 5) for n = 1 in two different ways. 

(i) The coefficients c^^ and can be calculated from Eqs. ([Ol ), ( |2r7|) in three dimen- 
sions according to 

F±(m) = -leTrZ^ei^ff (e±,47re±^Z„Z^-V3) + A(m,1), (3.8) 

where the Z-factors have the arguments Zi{u^ 1). The perturbative expression for f can 
be obtained for n = 1, 2, 3 from 

ff°)(e+,«o,3) = -^Z,-i(A) (3.9) 
4uo 
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where the renormahzation factor Z^IX) and its relation to the specific heat have been pre- 
sented in numerical form by Bervillier and Godreche |]28| and by Bagnuls and Bervillier 
29| , |36[] , see also Ref. ||T0[. For d = 3 their renormalized coupling A is related to our uq via 



the renormahzation factor Z-^^X) according to ^0^^+ = —^ttXZ^^X) as noted in Ref. |]T0 
The perturbative expression of f for n = 1 can be determined according to 

52 - 



^(2,0) 



(^_,Mo,3) 



r_(^_,Mo,3) 



dr'r, \ ^ d 



-r-(^-,uo,3) 



(3.10) 



where rQ(^_,Mo) is given by Eq. (3.8) of Ref. |Tl|]. The Helmholtz free energy f„(^_,'U0)3) 
is given in numerical form by Eq. (3.15) of Ref. |Tl| where our f_(^_,Mo,3) is denoted by 
r_o(^-,Mo)- Our numerical results for Cp^ up to nine digits are presented in Table fT^. 
(ii) Alternatively the coefficients be determined via the relation 

SA^^P±{u)ft''\u) = (1 - 2C.(n))F±(n) + AB{u) - (3u{u)dF±{u) / du (3.11) 

as done previously ||TU|,|rT|. For the definition of P± and f^''^'' and for a derivation of Eq. 
( ^.11| ) we refer to Refs. and [Q. In the present context we need the contributions to P± 
and f±'^^ only up to O(-u^) as given in Table 4 of Ref. |TD] and Table 3 of Ref. |]rT| since B{u) 
is known only up to 0{u^) as well. (We recall that the coefficients of P_ are determined by 
those of P+ according to P_(m) = — |{1 + 2[1 — P+{u)\ — \C,r{u)} 0.) This calculation via 
Eq. ( p.ll| ) yields coefficients that agree with those obtained via Eq. ( p.8|) up to eight 
digits for c~p^ and up to seven digits for c^^. The slight differences between the results of 
the calculations (i) and (ii) are due to the fact that Z^{X), Z^{X) and r_ are available up 
to five-loop order only in numerical form. We consider the calculation (i) via Eq. ( |3.8| ) as 
slightly more reliable since fewer numerical operations are required than in the calculation 
(n) using Eq. ( p.ll| ). 

Since here we have used the perturbative contributions of B{u) up to five-loop order, 
the resulting higher-order coefficients given in Table ^ differ from those determined 



previously (see Table 4 of Ref. [0 and Table 3 of Ref. |TT[) where the approximation 
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B{u) = n/2 + O^u^) was used. Only our low-order coefficients CpQ, Cp^, Cp^, Cpi and Cp2 
agree with the previous ones [|1^,|11|]. The coefficients Cp^ with m > 1 differ considerably 
from the previous ones whereas the coefficients Cpj^ with m > 2 differ only by 0.2% (m = 3), 
0.1% (m = 4), and 2% (m = 5). 

Very recently the coefficients and Cpj^ have been determined analytically for general 
n up to three-loop order |^. The corresponding three-loop values Cp2 for n = 1, 2, 3 as 



well as Cp^ for n = 1, 2, 3 are taken from Ref. [0] and are included in Table |ITl| up to 
nine digits. No results for Cp^ were available in the previous literature for n > 1. The new 
information on Cp^ for n > 1 enables us to perform the ffist Borel resummations for F^{u) 



and F_(u) — F<(u) for n = 2 and 3 in Subsection HID below. 



C. Borel resummation of five-loop results for n = 1 

In order to study the effect of the new higher-order terms we have performed Borel 
resummations of the series for uF^{u) and for u[F^{u) — at the fixed point u* for 

the case n = 1. The method employed is the same as for B[u) in Sect. II. The parameter 
ranges turn out to be 

1.6 < a < 1.7 , 7.48 <h < 8.70 (3.12) 

for u*F_{u*), and 

1.4<a<1.7, 6.0 <6< 11.7 (3.13) 

for m*[F_(m*) -F+(m*)]. 

We have found that our present method does not yield a reliable estimate of the param- 
eters a and b for F+{u*) separately; this may be related to the fact that, unlike Cp^^, the 
coefficients (see Table [1II| ) do not have alternating signs for m < 3 (this alternation 
is predicted for the asymptotic large-order behavior P,|55|). F^{u) will be further studied 



elsewhere. In the application to amplitude ratios given below we shall not need F+{u* 
separately. 
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The resummation results are 

u*F^{u*) = 0.3687 ± 0.0040 (3.14) 

and 

- F+{u*)] = 0.4170 ± 0.0036 . (3.15) 



The previous approximate resummation results |T3,11| for n = 1 were u*F_[u*) = 0.3648 and 



u*[F_{u*) — F^{u*)] = 0.4170 with an error bar of about 1%. Thus our resummation results 
differ from the previous ones only by about 1% for F_ and by less than 0.1% for F_ — F+, 
confirming previous conjectures |iy,ll]. The parameter dp in the effective representation 



111 F_{u) = {2u) ^ — 4(1 + dpu) now becomes dp = —4.64 (compared to —4.04 in Ref. 



TT|). This leaves the solid line in Fig. 4 of Ref. [11] essentially unchanged. 



D. Borel resummation of three- loop results for n = 2 and 3 

While the previous two- loop result for F^{u) did not yet provide sufficient information 
for a controlled resummation procedure for n > 1 and thus did not yet lead to an error 
estimate, the new three-loop coefficients Cp^ have significantly improved the situation. On 
the basis of these three-loop results we have performed Borel resummations of the series for 
uF_{u) and for u[F_{u) — F^{u)] for the cases n = 2 and 3. The method employed is the 
same as in Section |T|. The parameter ranges turn out to be 

- 1.70 <«< 1.10, 1.5<6<4.3, n = 2, (3.16) 

- 1.90 <«< 1.15, 1.5<6<4.7, n = 3 (3.17) 



for u*F_(u*), and 



-0.80 < a < 0.95, 0.65 < 6 < 8.3 , n = 2, (3.18) 
-0.30 < a < 0.55, 1.10 <6< 7.5, n = 3 (3.19) 



for u*[F^{u*) — F^{u*)]. The resummation results are 
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(m*) = 0.384 ± 0.025 , n = 2 



(3.20) 



u*F_ (u*) = 0.387 ± 0.026 , n = 3 



(3.21) 



and 



n*[F_(M*)-F+(n*)] = 0.461 ±0.019, n 



2, 



(3.22) 



m*[F_(m*) -F+(n*)] = 0.498 ±0.011, n 



3. 



(3.23) 



In order to substantiate the reliability of these resummations of three-loop results we 
have repeated our resummation for n = 1 but now on the basis of three-loop results. In 
Fig. 4a we show the Borel resummation value for uF^{u) at u = 0.04 for n = 1 together 
with the corresponding error bar. For comparison, the smaller error bar resulting from the 
Borel summation on the basis of the five-loop results is also shown. We see that the central 
value of the three-loop Borel resummation is not far from the central value of the five-loop 
resummation (arrow in Fig. 4a). Analogous results are shown in Fig. 5a for u[F_{u) —F^{u)] 
at u = 0.04 for n = 1. These results provide confidence in our resummation procedure based 
on three-loop results. In particular they demonstrate that Borel resummations of three-loop 
results yield reliable results with error bars that significantly reduce the uncertainties of one- 
and two-loop calculations. Since the n-dependence beyond three-loop order is expected to be 
smooth and weak in the range 1 < n < 3 there is no reason to expect that the reliability of 
our resummation procedure is significantly different for n = 2 and 3. In Figs. 4b and 5b the 
corresponding resummation results (with error bars) on the basis of three-loop perturbation 
series are shown for = 2 at the same value u = 0.04 as for = 1 in Figs. 4a and 5a. They 
clearly demonstrate the significant improvement over the previous situation at two-loop 
order P,p!^ where no error bar could be determined in a convincing fashion. 



IV. APPLICATIONS 



We apply our results to the specific heat in the asymptotic critical region where it can 
be represented as 
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a 



1 + a^|t|^ 



B 



(4.1) 



with the Wegner exponent A. We consider the universal amphtude ratios [|l| /A~ 
and la~. The former can be expressed in terms of B* = B{u*) and = F±{u*) in three 
dimensions as M 



A- U- 



1 — a 



F*-F* 
AuB* + aF* 



(4.2) 



where 



2uPl 



(4.3) 



6_ (3/2) - 2uP^ 

with = P^{u*). This expression for A'^ /A~' has several independent sources of inaccura- 
cies: (i) The values for the critical exponents a and z/, (ii) the values for P^, (iii) the values 
for 5*, (iv) the values for — F^ and F*. In evaluating this expression for A^ /A~ for 
= 1, 2, 3 we shall take (i) the most reliable values for a and v that are presently available, 
(ii) the Borel resummed values for P| based on five-loop results for n = 1, 2, 3, (iii) the 
Borel resummed values for 5* for n = 1, 2, 3 as given by Eqs. ( p.34|) -( ^.36| ), and (iv) the 
Borel resummed values for F^ as given in Eqs. ( p.l4|) , ( p.l5| ), and ( p.20| )-( ^.23| ). For the 
critical exponents a and u we take 



a = 0.1070 ±0.0045, i/ = 0.6310 ± 0.0015 , n = l, 

a = -0.01285 ±0.00038, z/ = 0.67095 ± 0.00013 , n = 2. 

a = -0.1150 ±0.0090, z/ = 0.7050 ± 0.0030 , n = 3. 



(4.4) 
(4.5) 
(4.6) 



according to Refs. |T^, and |3^, respectively. Borel summation values for have been 
calculated previously [110]. These calculations, however, were based on five-loop coefficients 



that were derived from the five-loop results of Refs. [1/] and |T8| which were corrected 



by Ref. |[T^. Taking into account the latter corrections we have derived corrected five-loop 
coefficients cps of the power series P+{u) = X^cT cpmu"^- Our new values for cps are (compare 
Table 4 of Ref. M) 



Cp5 = -1026631.56, n = l, (4.7) 

Cp5 = -1763840.20 , n = 2 , (4.8) 

Cp5 = -2767879.03 , n = 3 . (4.9) 

For the corresponding Borel resummed values of P| we have obtained 

= 0.7568 ± 0.0044 , n = l, (4.10) 

= 0.7091 ± 0.0045 , n = 2, (4.11) 

P* = 0.6709 ± 0.0039 , n = 3. (4.12) 



The values obtained from Eq. (3.2) and Table 1 of Ref. are by about 0.7% smaller 
than the values in Eqs. (|4:T0|)- (|4l^) . 

If we use the central values of the critical exponents given by Eqs. ( [4.4| )-( ^l6D and collect 
the results of Eqs. (CT)-(^^. ( CT) , ( pTTSD , (p:20D-(^:23D, and (CT)-(ra) 



we arrive at the values for A'^/A" according to Eqs. ( f4.2| ) and ( [4.3| ) as given in Table fV\ for 
n = 1, 2, 3. For comparison this Table also contains the results of previous calculations as 
well as experimental or numerical results. We see that there is good agreement between the 
previous and our results for n = 1 and n = 3. For n = 2, however, our result is significantly 
more accurate and agrees well with the high-precision experimental result for ^He near the 
superfiuid transition obtained from a recent experiment in space 0. 

The smallness of our error bar for A'^/A" for n = 2 is due to the small value of a for 
n = 2 which, according to the structure of Eq. ( [4.2|) , suppresses the error of F^. Exploiting 
this structure, i.e., separating exponents from amplitude functions, is a particular advantage 
of our d = 3 formulation of field theory. This structure was not explicitly taken into account 
in the previous e-expansion analysis for n = 2 The previous result obtained from the e 
expansion up to O(e^) does not agree with the experimental result for n = 2 within the 
error bars. The results 1.05 and 1.58 of Ref. |0 for n = 2 and n = 3 are partly based on the 
one-loop form of F_ {u) for which no error bar was available in Ref. . 

The expression of a^/cic is more complicated and depends on the derivatives of the 
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functions F^{u), F^{u) — F+{u), P+{u), B{u), and C-i^) at the fixed point, as specified in 
Eqs. (4.24), (5.19) and (5.21)-(5.24) of Ref . 0. We have performed Borel resummations for 
these quantities on the basis of our new five-loop resuhs for n = 1. For /a^ we obtain 
the value 1.0 ±0.1. This is consistent with the previous result ^0[ 0.96 ± 0.25 for n = 1. 
Since for n > 1 only three-loop results are available (which would still yield larger error bars 
for the derivatives of the various functions mentioned above) we postpone corresponding 
calculations of /a~ for > 1 until four-loop results become available. 

Finally we note that not only the fixed point values B* and but also the entire 
functions B{u) and F±{u) are of physical relevance. They are needed (i) for determining the 
effective renormalized coupling from experimental data of the specific heat , and (ii) for 
representing the specific heat in the entire nonasymptotic critical region well away from Tc 
according to Eqs. ( |3.3| )- (|3r5| ). We also note that -B* and F^ enter various other important 
universal ratios, e.g., related to the superfluid density ps of ''He 0. For this ratio we 
refer to a recent paper where the three-loop amplitude functions of the specific heat 
below Tc and of the order parameter are derived for general n. 

After completion of the present work we learned of the preprint [hep-ph / 9 71 0346| by B. 
Kastening, "Five-Loop Vacuum Energy Beta Function in 0^ Theory with 0(n)-Symmetric 
and Cubic Interactions" where the perturbative terms of a function equivalent to B{u) have 
been calculated up to five- loop order. These terms agree with ours in Eq. (|2.21| ). We have 
been informed by the author that this paper is accepted for publication in Phys. Rev. D. 
At the end of this work it is asserted that in our paper we use the e expansion. Our work 
does not use the e expansion to calculate the amplitude functions F±{u) and the amplitude 
ratios /A" and /a~ . Only the fixed point value u* has been determined via a Borel 
resummation of the e expansion series for u*. 
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APPENDIX A: DERIVATION OF THE ADDITIVE RENORMALIZATION A{u, e) 



The following expressions 



give the ultraviolet d = 4 poles of the diagrams 



shown in Fig. ^defined by K R' {d'^ I^'^^ / dr^) according to the standard notations, see e.g. Ref. 
|TE|. Here R' denotes the incomplete ultraviolet i?-operation which subtracts sub divergences, 
and K denotes the operation of taking pole parts. We use subscripts (a.b) for the pole terms 
on the right-hand-sides of (|AT|)- ([A8|) that are identical with the numbers associated with the 
three- and four-loop diagrams of Ref. |jT6| (the first number a in the brackets indicates the 
number of loops and the second number b indicates the consecutive number of a diagram in 
Table 1 of Ref. [|16|). The subscripts on the right-hand sides of Eqs. ( |A9|) -( [AT^ ) correspond 



to the numbers of the five- loop diagrams of Ref. . We have multiplied the left-hand sides 
by factors (IGvr^)"* in m-loop order because of the definition of the bare four-point coupling 
Wn'^ go/24: in Refs. 

The one- and two-loop expressions read 

1 



(o2 



-- J ln(ro 



{IGnyKR' 
the three-loop expressions read 

/ 02 



drn 
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(2.2) 



{16nr KR' /(^)(ro,6) 
{16nrKR' /(^)(ro,6) 



21, 



(3.2) 



J(3) 



8/(3.4) + 12/(3.9) = J^'^ 



the four-loop expressions read 

/ 92 



(4.5) = J^'^ 



aG^rKR'l^^^I^'\ro,e)j =2 

(16vr^r/r/?'(|^/(^)(ro,e)) =0 

{WtcYKR' (^I^ /(''(ro, e)^ = 4/(4.9) + 6/(4.11) = J^'^ 

{lQ7r')^KR' (^1^ /(8)(ro, e)^ = 24/(4.12) + 6/(4.18) + 12/, 



(4.19) 



J(8) 



(Al) 
(A2) 



(A3) 
(A4) 



(A5) 
(A6) 
(A7) 
(AS) 
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and the five-loop expressions read 

/ Q2 



/(9)( 



il6.r KR' 

il6.r KR' /(-)(ro,e 
il6nr KR' /(-)(ro,e; 
ilGny KR' (^I^''\ro,e: 



(16.^)^i^i?'(^/(-)(ro,e) 



(167r2firi?'(^/(i«)(ro,e) 



(9ro 
(9ro 



2/nr = J^'^) , 
= , 
= , 
2/121 = , 
4/7 + 6/120 = J^^^^ 
2/8 = , 



(A9) 
(AlO) 
(All) 
(A12) 
(A13) 
(A14) 



2/14 + 12/15 + 12/19 + 24/21 + 4/23 + I8/106 = J^''^ , (A15) 
8/79 + 16/88 + 32/95 + I6/100 = J^''^ , (A16) 

2/81 + 8/84 + 4/99 = /^^), (A17) 

8/32 + 32/47 + 4/73 + 8/93 + 8/98 + 4/109 + 8/116 = J^''^ . 

(A18) 



The pole terms up to four loops are ||T6| 

/ 



'(1.1) 



'(2.2) 



'(3.2) 



2 

e 



2_ _ 3_ 

3e2 ~ 4^ 



/(3.4) = — - — 



'(3.9) 



'(4.5) 



'(4.9) 



+ 



3e3 3e2 3e ' 
16 

4 3 

~3^^ 2^' 
16 16 4 



/(4.11) = ——r + -r^ 



3e4 + 3e3 3e^ 



(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
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'(4.12) 



'(4.18) 



4 10 

3?^ 3? 



13 11-6C(3) 
3e2 6e 



8 8 4 2C(3) 

\ \ h 

3e4 ^ 3e3 e 

1 7 



2 

3^3 



-^(4.19) = -— + — - 



12e 



8 

(J 






3^" 
8 


3 




3^ " 
4 


3 


5 


15e3 
2 


10e2 
13 


~ 96e ' 
857 


15e3 
8 


20e2 
7 


^ 960e ' 
25 215 


15e4 
16 


7 


6e2 96e ' 
11 157 

60e2 ' 320e' 
5 


15e4 
4 


~ 5^ " 
3 


15e3 


10e2 


~ 96e ' 



48C(3) 48C(3) + 24C(4) ^ 14C(3) + 9C(4) - 16C(5) 



5e3 



5e2 



and the five-loop pole terms are [18 

h 
h 

/l4 
/l5 
/l9 
/2I 
I23 
-^32 
-^47 
-^73 
-^79 
/81 
-^84 
-^88 
-^93 
-^95 
-^98 
-^99 

hoo 



5e 



15e5 5e4 
16 8 



12 ^ _6 ^ 18C(3)-45 ^ 146 - 90C(3) - 9C(4) 



15e5 3e4 
16 16 



28 



5e2 ■ 30e 
6 + 4C(3) 32 - 12C(3) - 18C(4) 



5e2 



30e 



8 4 + 4C(3) 7 + 6C(3)-12C(4) 

5e5 5e4 5e^ 5e^ 5e 



16 16 

3? ~ 3? 
20 



8 4 - 4C(3) 



3e3 
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3? ~ 3^ ^ 3e3 ~ 
16 8 28 



44 - 24C(3) 
12^2 ' 

6 - 12C(3) 16 - 18C(4) 



15e5 3e4 
8 52 

5^ ~ " 
16 



15e4 



+ 



15e3 5e2 15e 

34 116 -24C(3) 56 - 44C(3) + 6C(4) 

15^ ^ 60? We 

28 6 - 12C(3) 16 - 18C(4) 



15e5 3e4 15e3 

4 14 19 

+ 



15e4 15e3 15e2 
4 2 7 

3^~^^ 6^' 
4 _ 6 1 81-48C(3) 
5e4 ~ 5^ ^ 5e2 ^ 16(k 



5e2 15e 
386 + 768C(3) 

960e ' 



(A27) 
(A28) 
(A29) 

(A30) 
(A31) 
(A32) 
(A33) 
(A34) 
(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 
(A41) 
(A42) 
(A43) 
(A44) 
(A45) 
(A46) 
(A47) 
(A48) 



24 



-^106 
Il2l 



64 _ 32 
15e5 ~ 5e4 
16 16 

5^ 
4 

3^ 



15e5 



+ 



15e4 
32 

32 



32 



+ 



5e3 
16 

5^ 
32 

15^ 



3e3 



16 2 

15^ ~ 15e ' 

52- 108C(3) 202 - 168C(3) - 18C(4) 
15^2 

250 - 96C(3) 

120e ' 



3e5 3e4 
32 32 8 

3?~ 3?^3^ 



(A49) 
(A50) 
(A51) 
(A52) 
(A53) 
(A54) 



In Eqs. ( |A37|) and ( |A5 1| ) the corrections found in Ref. ||T9[ have been taken into account. 
Note that in Eqs. ([A19| )-( |A54|) we have used e = 4 — whereas in Refs. |[T6| . |18| e denotes 
(4-rf)/2. 

Eqs. ([A1|)-( |A5^ ) determine the additive renormahzation according to Eq. ( |2.7| ) as 



A{u,t) 



n 8n(n + 2) 



i=3 



i=5 



18 
i=9 



(A55) 



The overall factor of 2 in Eq. ( [A55| ) arises from the d = 4 value of the factor (Ai^/A) ^/16 



TT 



which is needed to obtain A{u,e) from Eqs. ( |A1|) -( |A5^) according to Eq. ( p.TI) . The renor- 
malized coupling u enters Eq. (|A55|) in the form u/2; the factor 1/2 takes into account that, 
near d = A, uq = Afu + 0{u^) = Sn^u + 0{u^) = HIQtt^u + 0{u^)] [see Eq. (U)]. 
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APPENDIX B: Z FACTORS 



In deriving the coefficients of the perturbation series of the quantities -F±('u), P±('u), 
and /±'°''('u) we needed the Z-factors Z.^, Z^, and calculated previously [p!7|-[T9[| up to 
five-loop order. Since their explicit form is not available in the literature we present them 
here explicitly. They read 

5 



Z,[u, e) = 1 + 5] Zf^\e)u'^ + O(m^) 

m=l 
m=l 



(Bl) 
(B2) 
(B3) 



m=l 



with the following coefficients in m-loop order: 
Coefficients of ZJu, e): 





4(n + 2) 
e 


Z?\e) 


= 4(n + 2) 




16 , 
= y (^ + 2) 




2 , 

= — (n + 2 

3 ^ 



4(n + 5) 



15n + 111 



-278 - 61n Vlr? + 132n + 360' 



(B4) 
(B5) 

(B6) 



288C(4)(5n + 22) + 48C(3)(3n2 + lOn + 68) + 31060 - + 7578n 



1152C(3)(22 + 5n) + 1236^2 + 23580n + 74616 ^ 16(245^2 + 2498n + 6284) 



-192 



(n + 5)(2n + 13)(n + 6) 



(B7) 



Zf\^) 



15 



:^+2) 



(9600C(6)(55n + + 186) + 768C(5)(-5n2 + 72 + 14n) 
+ 288C(4)(29n2 + 2668 - 3^^ + 816n) + 768C(3)2(-2n2 - 145n - 582) 
+ 48C(3)(8208n + 17n^ + 940n2 + 31848) + 21n^ + 4525471^ + 1077120n 
+ 3166528) ^ - (30720C(5) {2n^ + 186 + 55n) + 576C(4) (5n + 22) {n - 22) 
+ 96C(3)(27n3 + Vll^n^ + 14456n + 45448) - 9871^ + 277280^^ + 3073376n 
+ 7449712)^ + (2304C(3)(13n + 74) (5n + 22) + 21576^^ + 685192n2 
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^o. ^ -. ^ . .r.r..r.\ 1 32(307976 + 3175271^ + 172176n + 1933n=^ 
+ 5017312n + 10459360)— ^ . 



+ 



384(5n + 34) (n + 6) (n + 5) (2n + 13) 



Coefficients of Z^{u, e): 



7(1 
7(2: 
7(3 

7(5 



0, 



4(n + 2) 



= -(n + 2)(n + 8) 



3 

2(n + 2) 



1 4 



5(n2 - 18n - 100) 4(71^ + 234 + 53n) _ 16(n + 8)^ 



-^(» + 2) 



( - 1152C(4)(5n + 22) + 48C(3)(-6n^ + 184 + n^ + 64n) 
- 22752n - 39n^ - 296n^ - 77056)- 



+ 



2304C(3)(5n + 22) - 60n^ + 135488 + 42440n + 1844/1^ 



16(n + 8)(3n2 + 269n + 1210) ^ 192(n + 8)^ 



Coefficients of Zu{u, e): 



^ 16 



(n + SV 5n + 22 



+ 



96C(3) (5n + 22) + 942n + 2992 + 3571^ 16(n + 8) (17n + 76) 



24(n + 8)^ 



zi:^ 



16 

T 



+ 



(480C(5)(2n2 + 55n + 186) - 72C(4)(n + 8)(5n + 22) 

24C(3)(63n^ + 764n + 2332) + 20624n + 1640n^ - 5n^ + 49912)- 
480C(3)(5n + 22) (n + 8) + 67424n + 153088 + 7736/1^ + 172n^ 

16(55n + 248)(n + 8)2 48(n + 8)^ 



4 
15 



(6912C(7) (25774 + 9261n + 686n2) - 28800C(6)(n + 8)(2n2 + 55n + 186) 
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+ 768C(5) (165084 + 7466n^ + 305n^ + 66986n) - 288C(4) (62656 + 4084n^ + 28084n 

+ 189n^) + 2304C(3)^(3264 - 59n^ - 6n^ + 446n) + 48C(3)(1264n^ - 13n^ + 1312864 

+ 551032n + 67432n^) + 20429248n + 2518864n^ + 195n^ + 40148480 + 39230n^) - 

- (99840C(5)(n + S){2n^ + 55n + 186) - 14976C(4)(5n + 22) (n + 8)^ 

+ 3456C(3)(91n^ + 15436n + 34144 + 2196n^) + 63219712n - 800n^ 

+ 420800n^ + 117768192 + 981171271^) 4 

66048C(3)(5n + 22) (n + 8)^ + 32(n + 8)(733n3 + 4018671^ + 353392n + 803328) 

512(193n + 875)(n + 8)3 3840(n + 8)51 ,^ , 
-4 + ■ (B18) 
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FIGURES 

FIG. 1. Vacuum diagrams up to five-loop order determining the Helmholtz free energy Fq for 
Mq = 0, ro > 0. Diagrams (6), (10) and (11) do not contribute to A{u,e). The pole terms derived 
from the vacuum diagrams are given in (|AlD -( A54 ) of App. A. 



FIG. 2. Partial sums Bm{u) = E-m=o Cs^^*"" of B{u), Eq. < ^^ , as a function of u for 



n 



from M = 1 (two- loop order) to M = 4 (five-loop order). Also shown is the Borel resummed 
result (solid line) which deviates from the two-loop result i?i = 1 by only 0.5 % at the fixed point 
u* = 0.0362. 



FIG. 3. Borel resummation result for the function B{u), Eq. ( 2.21| ), for n = 2 (solid line) 



obtained by interpolation between the resummed values of B(u) at Uk = ku*/W, k = 1, . . . , 10 
of the renormalized coupling u in the range < n < = 0.0362, with error bars. Also shown 
is the three- loop result B2{u) = 1 -|- 24?i^ (dashed line), compare Fig. The two- loop result is 
Bi = 1. The Borel values B{uk) are 1.000233, 1.00073, 1.0013, 1.0019, 1.0026, 1.0032, 1.0037, 
1.0043, 1.0048, 1.0053 for = 1, . . . , 10, respectively. 

FIG. 4. Amplitude function F^{u), multiplied by the renormalized coupling u, in three dimen- 
sions as a function of u in 1-, 2- and 3-loop order (solid lines) for n = 1 (a) and n = 2 (b) and in 
4- and 5-loop order for n = 1 (a) (dashed lines). The dot with error bars at u = 0.04 indicates the 
result of the Borel resummation on the basis of the three- loop series for uF^ (u) . The small error 
bar (a) indicates the Borel resummation result on the basis of the five-loop series for uF-[u) at 
u = 0.04 for n = 1, the arrow indicates the corresponding central value at ti = 0.04. 
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FIG. 5. Amplitude function F-{u) — F+{u), multiplied by the renormalized coupling u, in three 
dimensions as a function of it in 1-, 2- and 3-loop order (solid lines) for n = 1 (a) and n = 2 (b) and 
in 4- and 5-loop order for n = 1 (a) (dashed lines). The dot with error bars at n = 0.04 indicates 
the result of the Borel resummation on the basis of the three-loop series for u[F^{u) — F+(u)]. The 
small error bar (a) indicates the Borel resummation result on the basis of the five-loop series for 
u[F-{u) — F+(ti)] at u = 0.04 for n = 1, the arrow indicates the corresponding central value at 
u = 0.04. 
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TABLES 

TABLE L Symmetry and group factors of the vacuum diagrams shown in Fig. ||. Diagrams 
(6), (10) and (11) do not contribute to A{u,e). 



loop order 


diagram (i) 


symmetry factor 5^*^ 


group factor G^^\n) 


1 loop 


(1) 


1 


n 


2 loops 


(2) 


3 


i(n2 + 2n) 


3 loops 


(3) 


36 


i(n^ + 4n^ + 4n) 




(4) 


12 


i(n2 + 2n) 


4 loops 


(5) 


432 


^{n'^ + 6n^ + I2n^ +8n) 




(7) 


576 


^{n^ + 4n2 + 4n) 




(8) 


288 


^(n^ + lOn^ + 16n) 


5 loops 


(9) 


5184 


^{n^ + 8n'^ + 24n^ + 32n^ + 16n) 




(12) 


10368 


^(n'^ + 6n3 + 12n2+8n) 




(13) 


6912 


^(n'^ + 6n3 + 12n2+8n) 




(14) 


6912 


^(n'^ + 6n3 + 12n2+8n) 




(15) 


2304 


i(n^ + 4n^ + 4n) 




(16) 


2592 


g^(n'' + 8n^ + 32n2 + 40n) 




(17) 


20736 


-^(n^ + 12n3 + 36n2 + 32n) 




(18) 


10368 


g^(5n3 + 32n2+44n) 
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TABLE II. Coefficients /j of the functions f3{u), Cr('u) and C(^(u) for i = 1,2, 3, respectively, 
and coefficients c^^ of B{u), compare Eqs. (2.20) and (2.21), for n = 0, 1, 2, 3 up to five-loop order 



(m = 6 for (3, m = 5 for d^r and C^, m = 4 for B). For compare Table 1 of Ref. 









Cr 




B 


n = 








8 












32 


-80 


-16 









-672 


3552 


128 









43989.9534 


-223152.607 


-8000 









-4166409.19 


18836823.8 


500639.112 









498653403.0 








n = 


1 





12 





1/2 






36 


-120 


-24 









-816 


6048 


216 


9 






56245.8519 


-413813.942 


-14040 


-761.422836 






-5632017.54 


37512804.7 


958294.321 


44244.7100 






708814936.0 








n = 


2 





16 





1 






40 


-160 


-32 









-960 


9024 


320 


24 






69029.7505 


-660870.017 


-21120 


-2256.06766 






-7268274.40 


63662497.1 


1566676.69 


141294.329 






956636505.0 








n = 


3 





20 





3/2 






44 


-200 


-40 






36 



-1104 12480 440 45 

82341.6490 -967074.371 -29000 -4653.13955 

-9075019.76 98265069.9 2333667.84 310944.846 
1243816220.0 
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TABLE III. Coefficients of F+(n) and F_(n) for n = 1, 2, 3 defined in Eqs. and 
(|3.7|) , respectively. For cj^„, m refers to corresponding to (m + l)-loop order whereas for Cp^, 
m refers to corresponding to m-loop order. The coefficients and Cp^ (three-loop order) 

are taken from Ref. [^] 





m 




^Fm 


n = 1 





-1 


1/2 




1 


-0 


-4 




2 


-22.6976284 


72 




3 


-722.742498 


-5189.75474 




4 


34775.5861 


433582.586 




5 




-47754702.5 


n = 2 





-2 


1/2 




1 


-16 


-4 




2 


-92.5270090 


64 




3 


-2430.86460 


-5918.07320 




4 


102469.659 




n = 3 





-3 


1/2 




1 


-30 


-4 




2 


-233.488142 


56 




3 


-5742.02976 


-6607.95641 




4 


204463.777 
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TABLE IV. Universal amplitude ratio /A of the specific heat, Eqs. 
2, 3 in three dimensions. 



in-M, for n = l, 



field theory 



present work 



previous work 



lattice series 
expansions 



experiment 



0.540±0.011 



0.524 ±0.010'' 



0.541 ±0.014*^ 



0.523 ±0.009' 



0.56-0.63^^, 0.53° 



1.056±0.004 



1.0294 ±0.0134'^ 
1.05^ 



1.08S 



1.054 ±0.001'" 
1.058 ±0.004' 
1.067 ±0.003j 
1.088 ±0.007^ 



1.51 ±0.04 



1.521 ±0.022^ 
1.58^ 



1.52S 



1.40' 



^Bervilher |32| 
^Bagnuls et al. 



^Liu and Fisher |39] 
'iBeysens et al. [M 



^Kumar et al. |41| 
^Schloms and Dohm 
^Hohenberg et al. 
^Lipa et al. |] 



'Lipa and Chui |43] 
^Singsaas and Ahlers |44] 



'^Takada and Watanabe |45| 
'Kornblit and Ahlers |4(:] 
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1 loop: 

(1) 



2 loops: 




(5) (6) (7) (8) 



5 loops: 




(9) (10) (11) (12) (13) 




(14) (15) (16) (17) (18) 



1 and 2 loop 
Borel 



5 loop 



n=2 



/ 



3 loop 



4 loop 







0.01 



0.02 
U 



0.03 



0.04 



